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ABSTRACT 

The collinear Lagrange points of the Sun-Earth sys- 
tem provide an ideal environment for highly sensi- 
tive space science missions. Consequently many new 
missions are planed by ESA and NASA that require 
satellites close to these points. For example, the 
SOHO spacecraft built by ESA is already installed 
in the first collinear Lagrange point. Neither uncon- 
trolled spacecraft nor escape motors will stay close to 
the Lagrange points for a long time. In case an oper- 
ational satellite explodes, the fragmentation process 
will take place close to the Lagrange point. Appar- 
ently a number of spacecraft will accumulate close 
to the Lagrange points over the next decades. We 
investigate the space debris hazard posed by these 
spacecraft if they explode and fall back to an Earth 
orbit. From our simulation we find that, as expected, 
about half of the fragments drift towards the Earth 
while the other half drifts away from it. Around 2% 
of the simulated fragments even impact the Earth 
within one year after the explosion. 



1. INTRODUCTION 

The hazards to manned as well as unmanned satel- 
lites by breakups of satellites and rocket upper stages 
in low Earth orbit (LEO) and geo-stationary or- 
bit (GEO) have been e xtensively discussed in the 
literature (jjchn, 199C|; iMcKnight fc Nagl, 1993 ; 
Fucke, 1993t |Jehn, 1995|; Matncy & Settecerri, 1997; 



Houchin et al. , 1997 y~Since the early 1970ies ( |Far 



^uhar, 1973) there are considerations to use another 
class of orbits: quasi-stable trajectories around the 
Lagrange points of the Sun-Earth system. Especially 
the collinear points Li and L2 (see figure |l|) are in- 
teresting fOT_sola^_ghj[sics and space science appli- 
cations (Farquhar, 199S). Since orbits around the 
collinear libration points are inherently instable, the 
dwell time of rocket bodies or uncontrolled satellites 
is short and thus breakups of these objects are less 
likely than in LEO or GEO. Breakups of malfunc- 
tioning operational satellites, however, create a cloud 
of fragments on the stable manifold. The stable man- 
ifold connects the locations of periodic moti on in the 
six-dimeii sional phase space of the satellite (Richard 
3on, 198C| ). 
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Figure 1. Illustration of the position of the La- 
grangian points of eguilihrium in the Sun-Earth sys- 
tem. The Sun is at the centre of the plane represent- 
ing the orbital plane of the Earth, and the horizontal 
axis is fixed to the Sun-Earth line. 



In order to model the fragmentation process, we ap- 
ply a fragmentation model that combines the frag- 
ment mass distribution by Bess (1975) with the 
Av distrib ution taken from the explosion model by 
Reynolds ( 1990| ). From these models, the differential 
fragment mass distribution of a low intensity explo- 
sion IS given by (see also Jehn ( |1990| )): 



ndm 



1.71 X 10-'*Mt exp(-0.02056V^) 
for m > 1936 g /-. x 

8.69 X lO-^Mt exp(-0.05756V^) 
for TO < 1936 g 



with Mt being the total mass of the satellite. As- 
suming an average fragment mass density of p = 
4.7 g cm^'^, the fragment diameter d is given by 



d 
2 



3to 
47rp 



(2) 



According to Reynolds, the Av of a fragment in 
[ms~^] is given as a function of fragment diameter d 
in fml : 



log Av = 



-0.0676 (logd)' 
-1.514 



0.804 log d (3) 



We simulate an explosion close to L2 using the mod- 
els described above. A total number of 820 frag- 
ments is generated, that drift off the point of equilib- 
rium depending on the direction and the magnitude 
of the Av imposed on the fragment by the explo- 
sion. Figure || shows the distribution of Av in the 
simulated fragment cloud. We assume that the frag- 
ments are distributed isotropically. After the sim- 
ulated breakup we propagate the fragments for one 
year, taking into account the gravity of the Earth, 
the Sun, the Moon, as well as the oblateness of the 
Earth and solar radiation pressure. In the follow- 
ing sections we describe the motion of the simulated 
fragments, as well as the potential threat they pose 
to operational satellites. 
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Figure 2. Histogram of the distribution of Av of the 
generated fragments (total number: 820). 



2. MOTION OF THE FRAGMENTS 

The dynamics of the fragments is dominated by 
the effective potential in the Earth-fixed frame of 
the restricted three-body problem. Thus, the frag- 
ment motion depends on the Jacobian constant that 
is determined by the breakup process. Before the 
breakup, the satellite is assumed to have the Ja- 
cobi constant of the L2 point of Co — —3.0008926 
in canonical units (distance unit = 1 AU, mass unit 
— Mq + , time unit such that period of Earth's 



orbit equals 2-k). For this value of C, the satellite is 
at the intersection of zero- velocity curves (ZVC), i.e., 
in equilibrium motion. Since the satellite is in rest 
with respect to the Earth-fixed frame, the breakup 
always increases the value of C. Because we consider 
lunar gravity, the Earth's oblateness, and radiation 
pressure, the Jacobi "constant" is not exactly con- 
stant along the fragment's trajectories. The value of 
C can be decreased by lunar perturbations (e.g. by 
close fly-bys) , or by the effect of the Earth's oblate- 
ness during close Earth encounters. Figure || shows 
the distribution of C for all 820 fragments 30 days 
after the breakup. 
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Figure 3. Histogram of the distribution of the Jacobi 
constant of the generated fragments (total number: 
820j. The dotted line indicates the initial value of 
the constant before the breakup. 



The breakup creates a rather narrow distribution of 
the Jacobian constant around the equilibrium value 
Co . Fragments with C < Co are confined either to a 
region close to the Earth or to outside a ring around 
the Earth's orbit. This restriction can be seen in 
figure ^ that shows the ZVCs for fragments with dif- 
ferent values of C. For C > Co the fragments may 
move freely along the Sun-Earth line, but not along 
the Earth's orbit. In general two classes of orbits 
are possible; out-bound and in-bound. Fragments 
on out-bound orbits move away from the Earth so 
that they have larger heliocentric semi-major axes 
and smaller mean motions. Consequently they drift 
into the positive y-direction, which is opposite to the 
direction of motion of the Earth around the Sun (see 
figure^ (a)). These fragments, which constitute 44% 
of the simulated objects, continue on independent 
heliocentric orbits and therefore pose no threat to 
Earth satellites. The other class of orbits are in- 
bound to the Earth that potentially exhibit close 
lunar or Earth encounters. If no close encounters 
occur, the orbits around the Earth are very unsta- 
ble, as shown in figure ^ (b). Mostly lunar fly-bys 
decrease the orbital energy (and thus the value of C) 





Figure 5. Trajectories (solid lines) of an explosion fragments in the Earth's orbital plane. In panel (a) a region 
11 X 11 AU, and in panels (h)-(d) a region 4 x 4 AU are shown with the Earth at the origin. The dotted circle 
indicates the orbit of the Moon and the dashed line the Sun-Earth line. The Sun is located at 149 x 10^ km 
along the positive x-axis. The cross-markers along the trajectories indicate the fragment's position separated by 
20 days. Panel (a) shows an Out-bound directed trajectory, panel (b) a typical in-bound trajectory, and panels 
(c) and (d) in-bound trajectories that are captured into bound orbits around the Earth. 



of some fragments, so that they stay in bound orbits 
about the Earth for more than one year (see figures 
1(c) and (d)). 



3. POTENTIAL HAZARDS TO 
OPERATIONAL SATELLITES 

In our simulation 56% of the fragments moved to- 
wards the Earth along the Sun-Earth line. Whether 
or not they pose a threat to operational satellites de- 
pends on whether they reach LEO or GEO distances, 
i.e. on their perigee distance r^r. Figure ^ shows the 
distribution of for all 820 fragments, propagated 
over one year. The peak at 1.5 x 10^ km indicates 
the 353 (44%) fragments on out-bound orbits that 
never approach the Earth closer than the distance 
of L2. The 467 (56%) fragments on in-bound orbits 
almost always approach the Earth closer than the 
lunar orbit. Figure ^ shows a zoom of the fragment 
perigee distribution inside the lunar orbit. The max- 
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Figure 4- Zero velocity curves (ZVCs) in the Earth- 
fixed system for different values of the Jacobi con- 
stant. A plane 0.1 x 0.1 AU around the Earth is 
shown, the Sun (not on the diagram) is at the ori- 
gin. The contour lines (solid for C > Cq, dashed 
for C — Co, dotted for C < C'o) represent ZVCs for 
C = Co, Co ± 0.001, Co ± 0.0005, andCo ± 0.0001. 
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Figure 6. Histogram of the perigee distribution of 
the fragments (total number: 820) one year after the 
breakup. 



imum of the perigee distribution is at 2 x 10^ km, 
far outside the orbit of most operational sateUites. 
While less than 60 fragments approached the Earth 
closer than GEO, no fragment actually intersected 
the GEO ring. In addition, a crossing of fragments 
from L2 of the GEO ring is only possible around the 
equinoxes when the GEO ring intersects the Sun- 
Earth line. 

Due to their low initial angular momentum (or 
transversal velocity) with respect to the Earth, frag- 
ments can actually impact the Earth's surface. In 
our simulation 17 fragments reached the LEO envi- 
ronment {H < 2000km), and 13 (1.6%) impacted the 
Earth close to Earth escape velocity (see figure @). 
Despite the fact that the out-of plane motion of tne 
fragments stays small, they can impact the Earth at 
any latitude, because the Earth's diameter small is 
compared to the distance to L2 (see figure ^ . 



4. CONCLUSION 

The use of the collinear Sun-Earth Lagrange points 
for solar physics and space science application is ex- 
pected to increase. Therefore, a breakup of a con- 
trolled satellite due to a malfunction can not be ruled 
out. Since the dwell time of uncontrolled satellites 
and rocket bodies close to the libration points is 
short, they are unlikely to contribute to the haz- 
ard to operational satellites due to breakups at these 
points. If an explosion of a satellite occurs close to 
L2, about 50% of the fragments move in-bound to- 
wards the Earth. Collisions of the fragments with 
GEO satellites are very unlikely and are possible 



only if the breakup occurs close to the vernal or au- 
tumnal equinox. A low percentage (about 2%) of 
the fragments reaches the LEO environment and im- 
pacts the Earth's atmosphere with velocities close to 
the Earth's escape speed of 11 kms~^. These Earth 
impacts can happen at any geographic latitude and 
longitude. 
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Figure 7. Zoom of the perigee range from to 4 x 
10^ km in figure The dashed lines indicate the 
distance range of the GEO ring and the lunar orbit. 
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Figure 8. Zoom of theverigee range from to 3i?0 
(Earth radii) in figure [j- The dashed lines indicate 
the boundaries of the LEO region and the hashed area 
indicates the extend of the Earth. 
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Figure 9. Trajectory (solid line) of an impacting 
fragment in the Earth's vicinity. The x-direction is 
the Sun-Earth line and the z-direction is perpendic- 
ular to the Earth's orbital plane. At the bottom a 
segment of the fragment's trajectory from an earlier 
close encounter can be seen. The final segment of 
the trajectory approaches the Earth from the upper 
left corner. 



